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Lunar Electrical Conductivity, Permeability, and Temperature 
from Apollo Magnetometer Experiments 
P. Dyal*, C. W. Parkin**, and W. D. Daily*** 

Abstract — Magnetometers have been deployed at four Apollo sites on the moon to 
measure remanent and induced lunar magnetic fields. Measurements from this 
network of instruments have been used to calculate the electrical conductivity, 
temperature, magnetic permeability, and iron abundance of the lunar interior. 
Global lunar fields due to eddy currents, induced in the lunar interior by magnetic 
transients, have been analyzed to calculate an electrical conductivity profile for the 

moon. From nightside magnetometer data in the solar wind it has been found that 

-4 

deeper than 170 km into the moon the conductivity rises from 3x10 mhos/m to 
_2 

10 mhos/m at 1000 km depth. Recent analysis of data obtained in the geomagnetic 
tail, in regions free of complicating plasma effects, yields results which are 
consistent with nightside values. Conductivity profiles have been used to calculate 
the lunar temperature for an assumed lunar material of olivine. In the outer layer 
(--170 km thick) the temperature rises to 1100°C, after which it gradually increases 
with depth to 1500°C at a depth of — 1000 km. Simultaneous measurements by 
magnetometers on the lunar surface and in orbit around the moon are used to 
construct a whole-moon hysteresis curve, from which the global lunar magnetic 

permeability is determined to be p, = 1.012 + 0. 006. The corresponding global 

18 3 -4 

induced dipole moment is 2x10 gauss-cm for typical inducing fields of 10 

gauss in the lunar environment. Lunar free iron abundance corresponding to the 

global permeability is determined to be 2. 5 + 2. 0 wt%. Total iron abundance (sum 

of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed 

compositional models of the lunar interior. For a free iron/orthopyroxene lunar 
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composition the total iron content is calculated to be 12. 8+1.0 wt%; for a free iron/ 

olivine composition, total iron content is 5. 5 + 1. 2 wt %. Other lunar models with an 

iron core and with a shallow iron- rich layer are also discussed in light of the measured 

global lunar permeability. The measured lunar remanent fields range from 3 gammas 

(v) minimum at the Apollo 15 site to 327 v maximum at Apollo 16. Simultaneous magnetic 

field and solar plasma pressure measurements show that the remanent fields at the 

Apollo 12 and 16 sites interact with, and are compressed by, the solar wind. Remanent 

fields at Apollo 12 and 16 are increased 16vand 32v, respectively, by a solar plasma 

-7 2 

bulk pressure increase of 1. 5 x 10 dynes /cm . Velocities and thicknesses of the 
earth’s magnetopause and bow shock have been estimated from simultaneous magneto- 
meter measurements. Average speeds are determined to be about 50 km/sec for the 
magnetopause and 70 km/sec for the bow shock, although there are large variations 
in the measurement for any particular boundary crossing. Corresponding measured 
boundary thicknesses average to about 2300 km for the magnetopause and 1400 km for 
the bow shock. 

INTRODUCTION 

In this paper we discuss the experimental results and analyses from magnetometer 
experiments deployed by astronauts on the surface of the moon in the Apollo program. 
Measurements from this network of instruments at the Apollo 12, 14, 15, and 16 
landing sites (shown in Figure 1) have been used to calculate the magnetic permeability, ' Fig. 1 
electrical conductivity, and temperature of the lunar interior. The fossil remanent 
magnetic fields have been measured at the four landing sites and the magnetic field 
interaction with the solar wind plasma has been measured at three of the sites. 

Simultaneous magnetic field measurements by the Explorer 35 lunar orbiting 
satellite and the Apollo surface experiments have been analyzed to determine the 
time dependence of eddy-current induction in the moon by solar Reid transients. 

This analysis has yielded a continuous electrical conductivity profile from which an 
internal temperature profile is calculated for assumed lunar material compositions. 

Similarly, magnetic field measurements obtained when the moon is immersed in the steady 
geomagnetic tail are used to construct a global hysteresis curve, from which the magnetic 



permeability of the lunar interior is determined. The iron abundance of the moon 
is then calculated for different compositional models of the lunar interior. 

The remanent magnetic field has been measured at the Apollo 12, 14, 15, and 16 
landing sites. The two photographs in Figure 2, taken by the Apollo 16 astronauts , 
show the two magnetometers used to measure remanent and induced fields at the 
Apollo 16 Descartes site. The remanent magnetic field provides a record of the 
magnetic field environment that existed at the moon 3. 7 to 4. 2 billion years ago at 
the time the crustal material cooled below its Curie temperature. The remanent 
field is found to be compressed by the solar wind plasma, and the properties of the 
compression are used to study the scale size and topology of the field. Simultaneous 
measurements obtained with this network of magnetometers on the lunar surface are 
also used to determine the thickness and velocity of the earth's bow shock and magneto- 
sheath. The positions of these boundaries relative to the lunar orbit are shown in 
Figure 3. 

ELECTRICAL CONDUCTIVITY AND TEMPERATURE OF THE LUNAR 

INTERIOR 

Electrical conductivity and temperature of the moon have been calculated from 
global eddy current response to changes in the magnetic field external to the moon. 
When the moon is subjected to a change in the external field, an eddy current field 
is induced in the moon which opposes the change (see Figure 4). The induced field 
responds with a time dependence which is a function of the electrical conductivity 
distribution in the lunar interior. Simultaneous measurements of the transient 
driving field (by Explorer 35) and the lunar response field (by an Apollo surface 
magnetometer) allow calculation of the lunar conductivity. Since conductivity is 
related to temperature, a temperature profile can be calculated for an assumed 
compositional model of the lunar interior. 

When the moon is in the solar wind, lunar eddy current fields form an induced 
lunar magnetosphere which is distorted in a complex manner due to flow of solar 
wind plasma past the moon. The eddy current field is compressed on the dayside 
of the moon and is swept downstream and confined to the “cavity” on the lunar nightside. 

Because of the complexity, early analysis included a theory for transient response of a 
sphere in a vacuum in oi’der to model lunar response as measured on the lunar 
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nightside (Dyal et al. , 1970a; Dyal and Parkin, 1971a; Sill, 1972i; a harmonic 
theory of a sphere totally confined by the solar wind was used to model the response as 
measured on the lunar dayside (Sonett et al. , 1971, 1972; Kuckes, 1971; Sill, 1972). 

Both the transient and the harmonic techniques have subsequently been further 
developed. Transient analysis has evolved to include effects of cavity confinement 
on nightside tangential data and to introduce analysis of magnetic step transients 
measured on the lunar dayside (Dyal and Parkin, 1973; Dyal et al. , 1973). Harmonic 
analysis has been developed with the purposes of eventually developing a dynamic response 
theory for the case of asymmetric confinement (Schwartz and Schubert, 1973) and of 
extending the data analysis to lower frequencies and accounting for diamagnetic effects 
in the solar wind (Kuckes et al. , 1974). 

Recently time- dependent poloidal response of a sphere in a vacuum has been 
applied to data measured in the geomagnetic tail where plasma confinement effects 
are minized. The poloidal response analysis has been used to determine the electrical 
conductivity and temperature profiles of the lunar interior. 

Electrical conductivity analysis: moon in solar wind plasma 

1. Lunar nightside data analysis. The lunar electrical conductivity has been 

investigated by analysis of the lunar response to transients in the solar wind magnetic 

field. The response, measured by an Apollo magnetometer on the nightside of the moon, 

is theoretically approximated by the response of a conducting sphere in a vacuum. The 

theory was developed by extending the work of Smythe (1950) and Walt (1951) for a radially 

varying lunar conductivity profile (Dyal et al. , 1972a). The measured response 

(illustrated in Figure 5a) is the average poloidal field response (for the radial surface Fig. 5 

field component) to a normalized fast- ramp decrease in the external field. Error 

bars are standard deviations of the measured responses. Chosen for the analysis is 

a ramp input function which falls from unity to zero in 15 seconds, a time characterizing 

convection of a solar wind discontinuity past the Moon. (For a 400 km/sec sola* wind, 

this time is 10-20 sec, depending on the thickness of the discontinuity and the inclination 

of its normal to the solar wind velocity). 

For a family of conductivity profiles, all of which monotonically increase with depth 
in the moon, the theoretical response to a fast ramp is calculated and compared to the 
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measured response. A particular set of these conductivity profiles yield response 
functions which pass within all data error bars of Figure 5a. These profiles define 
the shaded region of Figure 8 and are all consistent with the nightside response data. 

2. Lunar dayside data analysis. Theoretical solutions for an eddy current field 
totally confined to a sphere of homogeneous conductivity are derived from Maxwell's 
equations in Dyal et al. (1973). Figure 5b shows averages of normalized rising step 
transients measured on the lunar dayside in response to increasing step transients 
in the free-streaming solar wind (error bars are standard deviations). The over- 
shoot maximum is amplified by a factor of 5 over the external input field step change, 
by solar wind dayside confinement of the surface tangential field components. The 
data are fit by a lunar conductivity model with a homogeneous core of radius R c = 0. 9R^ 
and conductivity o~10 ^ mhos/m. This result is consistent with the nightside conductivity 
profile illustrated in Figure 8 to depths allowed by the duration of the response data 
which is shown in Figure 5b. 

Electrical conductivity analysis: moon in the geomagnetic tail 

The theoretical models outlined so far have all assumed spherical symmetry to 
describe lunar eddy current response to changes in the external field. However, the 
nightside and dayside analyses have used data taken when the moon is immersed in the 
solar wind plasma with asymmetric confinement of the inducing fields. The shortcomings 
of using spherically symmetric approximations to describe the induced lunar magneto- 
sphere, which is actually asymmetrically confined, have been pointed out in the 
literature for both the nightside vacuum approximation (see e. g. , Schubert et al. , 1973a) 
and the dayside totally-confined approximation (see, e. g. Dyal and Parkin, 1973) . Three- 
dimensional, dynamic asymmetric confinement presents a difficult theoretical problem 
which has not been solved at the time of this writing. Previous theoretical 
approximations of the asymmetric problem have included a two-dimensional approximation 
(Reiszetal., 1972); three-dimensional static theory for a point-dipole source, with 
substantiating laboratory data (Dyal and Parkin, 1 973); a three-dimensional “quasi-static” 
approach (Schubert et al. , 1973b) ; and a three-dimensional dynamic theory for one 
particular orientation of variations in the external field (Schwartz and Schubert, 1973). 

In order to circumvent this problem of asymmetry, recent analysis has considered lunar 

eddy current response during times when the moon is in the geomagnetic tail where plasma 
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interaction effects encountered in the solar wind (asymmetric confinement, 
remanent field compression, plasma diamagnetism, etc.) are minimal. 

1 Po loidal response of a sphere in a vacuum: theory. To describe the response of 

lunar sphere to an arbitrary input field in the geomagnetic tail, we define the magnetic 

vector potential A such that VxA = B and v* = 0. We seek the response to an 

input A B b(t), where b(t) = 0 for t<0 and b(t) approaches unity as t -*■ ». (Since the 

governing equations are linear, the response to a more general input is readily found 

by superposition.) The direction of A B ^ is taken to be the axis of a spherical coordinate 

system (r, 9 , cp) . If the conductivity is spherically symmetric, the transient magnetic 

field response has no cp component, and hence A = Ae^and 3/3 cp= o. Under these 

conditions (and neglecting displacement currents) the laws of Faraday, Ampere, and 

Ohm combine to yield a diffusion equation for the magnetic potential (in MKS units): 

• • 

V 2 X (r, 6 j t) = |J a (r) (r, 6; t) 

We will show in a later section on magnetic permeability that we may take M* = 
everywhere. Then, for t^O, the magnetic field must be continuous at the surface, so 
that A and 3A/3 r must always be continuous at r = R , the radius of the sphere . We 
also have the boundary condition A (0 , t) = 0 and the initial condition A (r, 9, o) = 0 inside 
the moon. Outside of the moon, where cr = 0, 


A = ABj, (|) b(t) sin0 + Aa, f(t) sin9 


( 2 ) 


The first term on the right is a uniform magnetic field modulated by b(t) ; the second 
term is the (as yet unknown) external transient response, which must vanish as r and 


Note that at r = R , where R is normalized to unity, 
m m 


*= AB„ sice 


(3) 


and 

3 A „ /b(t) \ 

gp = AB e sin 6 ^-y- + 2f(t)J . 


( 4 ) 





Therefore, at r = R = 1, _ . • . 

m !? - : « + 1 K .i* » ( t)) . 

' ( 5 ) 

Since the magnetic field is continuous at r = R , * s a boundary condition for the 
interior problem. Letting G(r, t) = A/A B sin 9 and G (r, s) be the Laplace transform 

Jlr 

of G, equation (1) becomes 
for the interior. The boundary conditions 


§g) = s\i Q a (r)G 
are 


at r=R and 
‘ m 


v 


at r = 0. 


G cO 


(?) 

(8) 


For a given a(r) and b(t), this system is numerically integrated to obtain G(r,s) in 

the range O^r^R . The function G< r, s) is then numerically inverse Laplace trans- 
m 

formed to find the characteristic transient response function f(t) for the system. The 
resulting f(t) is compared to the measurements and reiterated with a different function 
a (r) until the error between the calculated f(t) and the measured f(t) is minimized. The 
final o(r) is not unique; rather a family of c(r) is generated with the constraint that f(t) 


match the experimental data. 


2. Conductivity results: geomagnetic tail data analysis. Figure 6 shows an example 
' of a magnetic transient measured in the geomagnetic, tail. The data components are 
expressed in a coordinate system which has its origin on the lunar surface at the Apollo 
12 magnetometer site. The x-component is directed radially outward from the lunar 
surface, while the y-and z-components are tangent to the surface, directed eastward 
and northward, respectively. The external (terrestrial) driving magnetic field is 
measured by Explorer 35, whereas the total response field is measured on the lunar 
surface by the Apollo 12 magnetometer. 
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To analyze this and similar events, the Explorer 35 input field components are 
fitted numerically by functions which are sums of ramp inputs chosen to fit the 
external field data. For a given conductivity model of the moon, the theoretical time 
series response to an external field data set is numerically calculated using the 
theory outlined in the previous section. This calculated response is then compared with 
the measured Apollo 12 time series data. Conductivity profiles have been constrained 
to be monotonically increasing with depth into the moon. 

Figure 7 shows an example of calculated response for the Explorer 35 x-axis 
(radial) input function of Figure 6, using the conductivity profile illustrated in Figure 
7. Superimposed is the actual response, which is the Apollo 12 x-component of 
Figure 6. This conductivity profile yields the best fit of seventy profiles which have 
been run to date, although it is certainly not unique. The profile also yields 
theoretical responses which fit well for the measured tangential components of Figure 6 
and the components of "twenty other geomagnetic tail transients which have been 
processed to date. The analysis presented here is still preliminary; in the future 
many more transient events will be processed to determine a range of conductivity 
profiles consistent with a large data set. 

In the previous sections various techniques which we have utilized to investigate 
lunar electrical conductivity were surveyed. Figure 8 shows a summary of 
conductivity profiles which have been calculated. 

Lunar temperature profiles from conductivity analyses. Once the electrical conductivity 
profile of the moon has been determined, an internal temperature distribution can be 
inferred for an assumed lunar material composition {see Rikitake, 1966). For cases 
where electrical conductivity is independent of pressure to a first approximation, the 
conductivity of materials can be expressed in terms of temperature T as follows: 


a=?E i exp (-a^RT) (S) 

where a. are the activation energies of impurity, intrinsic, and ionic modes, expressed 
In electron volts; are material-dependent constants; and k is Boltzmann's constant. 

It should be emphasized that the electrical conductivity o (a, E, T) is a strong function 


b 


Fig. 7 


Fig. 8 



of the material composition; therefore, uncertainties in knowledge of the exact 
composition of the sphere limits the accuracy of the internal- temperature cal- 
culation. - . 

Laboratory results relating conductivity to temperature for various minerals, 
which are good geochemical candidates for the lunar interior, have been presented 
by many investigators (e.g. , England et at. , 1968; Schwerer et al. , 1972; Olhoeft et 
al. , 1973; Duba et al. , 1972; Duba and Ringwood, 1973). The recent work by Duba 
et al. (1972) on the electrical conductivity of olivine has been used to convert the 
electrical conductivity profiles of Figure 8 to the temperature profiles shown in 
Figure 9. Fig. 9 

MAGNETIC PERMEABILITY AND IRON ABUNDANCE IN THE 

MOON 

Magnetic permeability and iron abundance of the moon are calculated by analysis 
of magnetization fields induced in the permeable material of the moon. When the 
moon is immersed in an external field it is magnetized; the induced magnetization is 
a function of the distribution of permeable material in the interior. Under the 
assumption that the permeable material in the moon is predominately free iron and 
iron-bearing minerals, the lunar iron abundance is calculated from the lunar 

i 

permeability for assumed compositional models of the interior. Since the amount of 
iron present in the lunar interior should be consistent with the measured global 
magnetic permeability, the permeability in effect places a constraint on the physical 
and chemical composition of the moon's interior. 

Global magnetic permeability 

Deployment of Apollo magnetometers on the lunar surface has allowed 
simultaneous measurements of the external inducing field (by Explorer 35) and the 
total response field at the lunar surface (by an Apollo magnetometer). During times 
when eddy current induction fields and interaction fields are negligible, the totai 
response field measured at the surface by an Apollo magnetometer can be expressed: 


! 7 H = H 4 AfTH 


( 10 ) 
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e H is the external magnetizing field and M. is the magnetization field induced 

e permeable lunar material (see Figure 10). The relative magnetic permea- Fig. 10 

3 

r is M-" 1 + 4 TTk, where k is magnetic susceptibility in emu/cm . Since the 
ar magnetization M is known to be below the Explorer 35 magnetometer reso- 
a (Behannon, 1968) , it is assumed in the dual magnetometer analysis that 
orer 35 measures H alone. 

For the two-layer lunar permeability model illustrated in Figure 10 (which 
be referred to later when iron abundance calculations are considered), the total 
at the outer surface of the sphere is expressed 


B - +2 G) x + Hy(l-G) y + H z (l-G) z 

(2iyfl)(p 1 -l) - X^(n-l)(2p 1 +l) 

. (2^l)d 1+ 2) - 2x 3 (n-i)(n x -i) 


( 11 ) 

( 12 ) 


i 1] = an ^ u ‘j are re ^ at ^ ve Permeability of the shell and core, respectively. 

permeability exterior to the sphere is p = 1, that of free space; \ = R /R ; 

o cm 

nd R are radius of the core and the moon, respectively. Equation (11) 
m 

esses the total surface field in a coordinate system which has its origin on the 
c surface at an Apollo magnetometer site: x is directed radially outward from 
unar surface, and y and z are tangential to the surface, directed eastward and 
lward, respectively. 

A plot of any component of equation (11) will result in a B - H hysteresis curve. 

ition (32) relates the slope of the hysteresis curve to the lunar permeability. The 

age whole-moon permeability p is calculated from the hysteresis-curve slope by 

ng p = p = p in equation (3) : 

1 2 


p - 1 

G = - — ~ 

p+ 2 


(IS) 


The hysteresis-curve method of permeability analysis was first employed by 
and Parkin (1971a) to calculate the %vhole-moon permeability result 1. 03 + 0. 13. 
3 then the error limits have been lowered by processing a larger number of 
litaneous data sets and using more rigid data selection criteria (e.g. , Parkin 


. , 1973). 
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In the most recent dual-magnetometer results {Parkin et al. , 1974a, b), a 
hysteresis curve has been constructed using 2703 data sets (see Figure 11). Since Fig. 11 

the external magnetizing field is so small (~T0 gammas), the familiar ”S" shape of 
hysteresis curve degenerates to a straight line (Ellwood, 1934). The data have been 
fitted by a least squares technique which yields the slope best estimate of 1. 008 + 

0. 004. Using this value with the radial (x) component of equation (11) and equation (13) , 
the whole-moon permeability has been calculated to be a = 1. 012 + 0, 006 (2 o error 
limits). Both extrema are greater than 1. 0, implying that the moon, as a whole, acts 
as a paramagnetic or weakly ferromagnetic sphere. This result has been used to 
calculate the iron abundance of the moon as discussed in the next section. 

Russell et al. (1974)have recently made permeability calculations using data from 
a single magnetometer, the Apollo 15 subsatellite magnetometer orbiting at an altitude 
about 100 km above the moon. The results to date indicate that the relative permeability 
of the entire spherical volume enclosed by the satellite orbit is below 1.0, implying 
that the layer between the moon and the satellite orbit is diamagnetic. Whether such 
a layer exists is uncertain at this time; further investigation is required using both 
magnetic and plasma data. 

Lunar iron abundance 

Iron abundance calculations have been presented by various authors, in theoretical 
treatments based on geochemical and geophysical properties calculated for bodies of 
planetary size (Urey, 1962; Reynolds and Summers, 1969; Urey and MacDonald, 1971) 

rt 

or on measured compositions of meteorites (Wanke et al. , 1973). Recently the global 

lunar permeability measurement, determined from magnetic field measurements, has 

been modeled by a homogeneous paramagnetic rock matrix (olivine and orthopyroxene 

models are used), in which free metallic iron is uniformly distributed. Pyroxenes and 

olivines have been reported to be major mineral components of the lunar surface fines 

and rock samples (Nagata et al. , 1971; Zussman, 1972; Weeks, 1972), with combined 

2 + 

iron present at the paramagnetic Fe^ ion. The ferromagnetic component of lunar 
samples is primarily metallic iron which is sometimes alloyed with small amounts of 
nickel and cobalt (Nagata et al. , 1972: Pierce et al. , 1971). This free iron is thought 
to be native to the moon (because of its low nickel content) rather than meteoritic in 
origin (Strangway et al. , 1973a) . Orthopyroxene and olivine models are consistent with 



jcbemical studies (Urey et al. , 1971; Wood et al. , 1970; Ringwood and 

11 

sene, 1970; Green et al. , 1971) and geophysical studies (Toksoz, 1974). 

Since the susceptibility of free iron changes several orders of magnitude at 

: iron Curie temperature (T c ), a two-layer model has been used, with the 

'-e-shell boundary R at the Curie isotherm (see Figure 10). For R >R , 

c c 

T . Therefore, for R >R any free iron is ferromagnetic while at greater 
c c 

)ths where T>T , the free iron is paramagnetic. The Curie isotherm location 
c 

determined from the thermal profile used for a particular model. Three thermal 

•dels have been used in the calculations. For model profile T the Curie isotherm 

spherically symmetric and located at R /R = 0. 9. Shell and core temperatures 

c rn. q 

3 600°C and 1400°C, respectively. For the model profile T g the shell is 500 C, 

i the core is 1300°C, while the Curie isotherm boundary is at R /R = 0, 85. 

mperatures are 300°C and 700°C for shell and core of model profile T^, which has 

/R = 0. 7. In the outer shell there are both ferromagnetic and paramagnetic 
m 

atributions to the total magnetic permeability - 1 + 4 nk^. The susceptibility of 

shell is k = k + k . where k, is "apparent' 1 ferromagnetic susceptibility and 
1 1c la la 

is paramagnetic susceptibility. The ferromagnetic component is metallic free 

n, assumed to be composed of multidomain, noninteracting grains: the paramagnetic 
2 + 

nponent is Fe~ combined in the orthopyroxene or olivine rock matrix. The measured 
Tomagnetic susceptibility of the shell material is an apparent value which differs from 
intrinsic ferromagnetic susceptibility of the iron because of self-demagnetization of the 
i grains and the volume fraction of iron in the shell. For R<R C the lunar materia! is 
amagnetic only, with susceptibility k 2 = k 2c + k 2a ; k 2c is the contribution of paramagnetic 
mically combined iron and k 2a is the apparent susceptibility of free paramagnetic iron 
>ve the Curie temperature. 

From the magnetic properties c? lunar compositional and thermal models, lunar 
m abundances have been calculated for the moon which are consistent with measured 
otal permeability (Parkin et al. , 1974a, b). The results are summarized in 

gure 12. The minimum total iron abundance consistent with the hysteresis curve can Fig. 
Iculated assuming the whole-moon permeability corresponds entirely to ferromagnetic 
on in the outer shell where the temperature is below the Curie point. For this case 
e bulk iron abundance is 0. 9 + 0. 5 wt. %. It is noted that the susceptibilities of both 
ivine and orthopyroxene are about an order of magnitude too small to account for the 
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measured permeability without some ferromagnetic material present. 

3. Considerations of an iron core and iron-rich layer 

The whole-moon permeability has also been used to investigate the magnetic effects 
of a hypothetical iron core in the moon. Density and moment of inertia measurements 
for the moon limit the size of such a core to less than 500 km in radius (Toksoz, 1974). 
If this hypothetical iron core were entirely paramagnetic and the surrounding core were 
orthopyroxene of average temperature 1100°C the global permeability would be 1.0003. 
This value is small compared to the measured permeability of 1. 012 + 0. 006, implying 
that if such a small paramagnetic iron core exists, its magnetization is masked by 
magnetic material lying nearer to the surface. Therefore the hysteresis measurements 
can neither confirm nor rule out the existence of a small iron core in the moon. 

An iron- rich layer in the moon has been considered by several investigators (e. g. , 
Wood et al. , 1970; Urey et al. , 1971; Gast and Giuli, 1972). It is possible that early 
melting and subsequent differention of the outer several hundred kilometers of the moon 
may have resulted in the formation of a high-density, iron-rich layer beneath a low- 
density, iron-depleted crust. Constraints have been placed on an iron- rich layer by 
Gast and Giuli (1972) using geochemical and geophysical data (for example, measure- 
ments of lunar moments of inertia). One set of their models consists of high-density 
layers between depths of 100 km and 300 km. At a depth of 100 km the allowed layer 
thickness is 12 km; the thickness increases with increasing depth, to 50 km at 300 km 
depth. Also presented are a set of layers at 500 km depth. By using exactly the same 
considerations as were used in the iron abundance calculations, calculations have been 
made of the whole-moon permeabilities which would be expected from lunar models 
with these iron- rich layers (Parkin et al. , 1974a, b). The calculations indicate that all 
iron rich layers allowed by geophysical constraints as outlined by Gast and Giuli, 
if wholly above the iron Curie temperature and therefore paramagnetic, would yield 
global permeabilities of about 1. 00006. As for the case of a small lunar iron co^e, the 
magnetization field of such paramagnetic layers would be masked by ferromagnetic 
materials elsewhere in the moon, and the hysteresis curve measurements camneither 
confirm nor rule out these layers. This conclusion would particularly apply to the 
Gast-Giuli layers at 500 km depths, which are almost certainly paramagnetic. If the 
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Iron-rich layers are below the Curie temperature and therefore ferromagnetic, they 
would yield measured global permeabilities of about 3. 5. This value is well above . 
the upper limit for the actual measured permeability of 1. 012 + 0. 006, and therefore 
the Gast-Giuli layers can be ruled out if they are cool enough to be ferromagnetic. It 
is important to note that the high-density layers discussed by Gast and Giuli (1972) can 
be thought of as limiting cases and that there are innumerable less dense and thinner 
layers which are allowed by geophysical, geochemical and magnetic constraints. 

LUNAR REMANENT MAGNETIC FIELDS AND INTERACTIONS 
WITH THE SOLAR WIND 

The permanent magnetic fields of the moon have been investigated using surface 
magnetometer measurements at four Apollo sites, orbital measurements from 
Explorer 35 and two Apollo subsatellite magnetometers, and natural remanent magne- 
tization measurements of returned lunar samples. Lunar remanent field measurements 
by surface magnetization measurements have been considered elsewhere (Hinners, 1971; 
Nagata et al. , 1972; Strangway et al. , 1973a; Fuller, 1974). 

The permanent lunar magnetic fields were first measured in situ by the Apollo 12 
lunar surface magnetometer (LSM) which was deployed on. the eastern edge of Oceanus 
> rocellarum. The permanent field magnitude was measured to be 38 + 3 gammas 
and the source of this field was determined to be local in extend (Dyal et al. , 1970b; 
Barnes et al. , 1971). A remanent field this large was generally unexpected even 
though relatively strong natural remanent magnetism had been discovered in the 
Apollo 11 samples, and the explanation of the field's origin yet remains a central 
problem in lunar magnetism. Subsequent to this measurement of an intrinsic lunar 
magnetic field, surface magnetometers have measured fields at the Apollo 14, 15 and 
16 sites. Fields of 103 + 5 and 43+6 gammas, at two sites located about a kilometer 
apart, were measured by the Apollo 14 Lunar portable magnetometer (LPM) at Fra 
Mauro. A steady field of 3. 4 + 2. 9 gammas was measured near Hadley Rille by the 
\pollo 15 LSM. At the Apollo 16 landing site both a portable and stationary magneto- 
neter wer6 deployed; magnetic fields ranging between 112 and 327 gammas were 
measured at five different locations over a total distance of 7. 1 kilometers at the 
Descartes landing site. The Apollo 16 fields are the largest lunar fields yet 
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measured. A schematic representation of these measured field vectors is shown 
in Figure 13. All the vectors have components pointing downward except the one Fig. 13 

at Site 5 near Stone Mountain, which points upward. This suggests, among other 
possibilities, that the material underlying Stone Mountain has undergone different 
geological processes than that underlying the Cayley Plains and North Ray Crater. 

In fact, Strangway et al. (1973b) have proposed the possibility that the light colored, 
relatively smooth Cayley formation is magnetized roughly vertically; the difference 
in the vertical component at site 5 is explained as an edge effect at the Cayley Plains — 

Stone Mountain boundary. A summary of all remanent lunar fields measured ty the 

magnetometer deployed on the surface is given in Table 1. 

Interaction of the solar wind with the remanent magnetic field has been measured 

at the Apollo 12 and 16 landing sites. -The solar plasma is directly measured at the 

Apollo 12 and 15 sites (Clay et al. , 1972) and simultaneous magnetic field and plasma 

data show a compression of the steady field as a function of the solar wind pressure 

at the Apollo 12 and 16 sites (Dyal et al. , 1973). The nature of the correlation between 

magnetic field and plasma bulk flow pressures is shown in Figure 14, which shows Fig. 14 

data (combined from several lunations) at the Apollo 12 and 16 LSM sites. The plasma 

bulk flow pressure and the magnetic pressure are related throughout the measurement 

range, and the magnitudes of magnetic pressure changes are in proportion to the 

unperturbed steady field magnitudes at each site. 

Information on the scale sizes of the permanently magnetized regions near Apollo 
landing sites is given bv gradient measurements of the lunar suriace magnetometeis, the 
spacing of vector measurements over the lunar surface, the known interaction 
properties of these remanent fields with the solar wind plasma, and limits imposed 
by satellite measurements. The field gradient in a plane parallel to the lunar surface 
is less than the instrument resolution of 0. 13 gamma/meter at the Apollo 12 and 15 sites. 

At Apollo 14 a. field difference of 60 gammas was measured at two sites located 1. 1 km 
apart. Gradient measurements and the absence of changes in the permanent field at 
the sites after lunar module ascent have demonstrated that the field sources are not 
magnetized artifacts. 

The scale size of the Apollo 12 remanent field has been calculated from local 
gradient and Explorer 35 measurements to be from 2 km to 200 km (Dyal et al., 1972b). 

For the Apollo 16 field, portable magnetometer measurements over the lunar roving 
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vehicle traverse have shown that the scale size for the field is greater than 5 km; 
the Apollo 16 subsatellite magnetometer has shown no anomalous field attributable 
to the Descartes area at orbital altitude, implying a surface field scale size upper 
limit of 100 km. Therefore, the Apollo 16 remanent field scale size is between 5 
and 100 km. 
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PROPERTIES OF THE EARTH'S MAGNETOPAUSE AND BOW SHOCK 

Using simultaneous data from a net work of magnetometers on the lunar surface 
and in orbit around the moon, the velocity and thickness of the geomagnetic magnetopause 
and bow shock have been measured at the lunar orbit {Daily et al. , 1974). Elapsed-time 
data for shock and magnetopause motions as measured by magnetometers separated up 
to 10^ km from each other, indicate that these boundaries are nearly always in motion 
and can have highly variable velocities. The magnetopause has an average speed of 
about 50 km/sec but measurements vary from less than 10 km/sec up to about 150 
km/sec. Similarly, the bow shock has an average speed of about 70 km/sec but again 
there is a large spread in measured values from less than 10 km/sec to about 200 km/sec. 
Based on the calculated boundary speeds and the signature of the boundary in the 
magnetometer data, the thickness of the bow shock and magnetopause have been estimated. 

I 

The average measured magnetopause thickness is about 2300 km; however, 
individual magnetopause boundaries range from 500 km to 5000 km in thickness. The 
average bow shock thickness is determined to be about 1400 km, with a spread in 
individual values ranging from 220 km to 3000 km. 
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SUMMARY 


1. Lunar electrical conductivity and temperature 

The electrical conductivity of the lunar interior has been investigated by 
analyzing the induction of global lunar fields by time varying extralunar (solar 
or terrestrial) magnetic fields. An upper limit on the unipolar induction field has 
been determined (Dyal and Parkin, 1971b) which shows that at least the outer 5 km 
of the lunar crust is a relatively poor electrical conductor (<10 8 mhos/m) compared 
to the underlying material. Past conductivity analyses have all used magnetometer data 
recorded at times when global eddy current fields were asymmetrically confined by the 
solar wind plasma, using two basic approaches: a time-dependent, transient-response 
technique and a frequency-dependent, Fourier-harmonic technique. Transient analysis 
using lunar nightside data yields a conductivity profile rising from about 3 x 10~^ 
mhos/meter at 170 km depth in the moon to about 10~ 2 mhos/meter at 1000 km depth. 

Recent conductivity analysis has considered lunar eddy current response during 
times when the moon is in the geomagnetic tail in order to avoid the analytical problems 
posed by asymmetric solar wind confinement of the induced lunar magnetosphere. 

Preliminary results show that the following conductivity profile, though not unique, is 
compatible with input and response data: the conductivity increases rapidly with depth 

— 0 A 

from 10 mhos/meter at the surface to 10 mhos/meter at 340 km depth, then less 
rapidly to 2 x 10 mhos/meter at 870 km depth. This conductivity profile is converted 
to a lunar temperature profile using the data of Duba et al. (1972) for olivine, and 
compared to the thermal profile results from other analyses. 

2. Lunar magnetic permeability and iron abundance 

Simultaneous measurements by lunar magnetometers on the surface of the moon 

and in orbit around the moon have been used to construct a whole-moon hysteresis curve, 

from which the global lunar relative magnetic permeability is determined to be 1. 012 + 0. 006. 

The global induced magnetization dipole moment corresponding to the permeability 

22 

measurement is 2x10 II (where H is magnetizing field in gauss). For typical geomagnetic 
tail fields of 10 gauss, the corresponding induced dipole moment is 2 x 10 18 gauss- 
cm 3 . Both error limits on the magnetic permeability value are greater than 1.0, implying 
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that the moon as a whole is paramagnetic and/or weakly ferromagnetic. Assuming 
that the ferromagnetic component is free metallic iron of multidomain, noninteracting 
grains, the free iron abundance in the moon is calculated to be 2. 5 + 2, 0 wt. %. Total 
iron abundance in the moon is determined by combining free iron and paramagnetic iron 
components for two assumed lunar compositional models. For an orthopyroxene moon 
of overall density 3.34 g/cm 3 with free iron dispensed uniformly thoughout the lunar 
interior, the total iron abundance is 12. 8 + 1. 0 wt. %. For a free iron/ol iv i ne moon 
the total iron abundance is 5. 5 + 1. 2 wt. %. 

3. Lunar remanent magnetic fields 

Direct measurements of remanent fields have been made at nine sites on the lunar 
surface: 38 v at Apollo 12 in Oceanus Procellarum;103 v and 43 vat two Apollo 14 sites 
separated by 1: 1 km in Fra Mauro; 3 v at the Apollo 15 Hadley Rille site; and 189 v , 

112v , 327 y, 113 y , and 235 y at five Apollo 16 sites in the Descartes region, over a 
distance of 7. 1 km. Simultaneous data from Apollo surface magnetometers and solar 
wind spectrometers show that the remanent fields at the Apollo 12 and 16 sites are 
compressed by the solar wind. In response to a solar wind dynamic pressure increase 
of 1. 5 x 10*" 7 dynes/cm 2 , the 38 gamma remanent field at the Apollo 12 LSM site is 
compressed to 54 gammas, whereas the field at the Apollo 16 LSM site correspondingly 
increases from 235 gammas to 265 gammas. Scale sizes of fields at the Apollo 12 and 16 
sites have been determined from properties of the remanent field-plasma interaction and 
orbiting magnetometer measurements. The Apollo 12 scale size L is in the range 
2 km£L £200 km, whereas for Apollo 16, 5 km < L < 100 km. 

Thei'e is no direct evidence that a global permanent field exists on the moon at 
present. Analysis of Apollo subsatellite magnetometer data (Russell et al. , 1973) has 
yielded an upper limit of 4. 4 x 10 13 gauss-cm 3 for the global permanent field dipole 
moment. 

Measurements by lunar orbiting satellite magnetometers and the Apollo surface 
magnetometers, and remanence in the returned samples, have yielded strong evidence 
that the lunar crustal material is magnetized over much of the lunar globe. The origin of 
the lunar remanent field remains an enigma. Possibilities are generally grouped under 
three classifications; a strong external (solar or terrestrial) field, an ancient intrinsic 
field of global scale, and smaller localized field sources. 
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4. Velocity and thickness of the earth’s magnetopause and bow shock 


Velocities and thicknesses of the earth's magnetopause and bow shock have been 
estimated from simultaneous magnetometer measurements. Average speeds are 
about 50 km/sec for the magnetopause and about 70 km/sec for the bow shock, with 
large spreads in individual measured values. Average thicknesses are about 2300 km 
for the magnetopause and 1400 km for the bow shock, again with large spreads in 
individual measured values. 
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XftTjle i. SnHuns.i'y of Lunar Surface Remanent Magnetic Field Measurements * 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Apollo magnetometer locations and measurements on the lunar surface. 
Maximum remanent magnetic fields measured at each landing site are 
shewn. The insert shows an orbit of the Explorer 35 spacecraft, projected 
onto the solar ecliptic plane. 

Apollo 16 magnetometers deployed on the moon at the Descartes landing site. 

(a) Lunar surface magnetometer (LSM). Sensors are at the top ends of the 

booms, approximately 75 cm above the lunar surface, (b) Lunar portable 

magnetometer (LPM), deployed during a magnetic field measurement by 

astronaut Young. Sensors are located on top of tripod at far right in photograph. 

Magnetic environment of the moon during a lunar orbit, with emphasis on 

the geomagnetic tail region. The plane of the lunar orbit very nearly 

coincides with the ecliptic plane of the earth 's orbit. The earth's permanent 

dipole field is swept back into a cylindrical region known as the geomagnetic 

-4 

tail; at the lunar distance the field magnitude is~10 gammas (10 gauss). 
Substructure of the tail includes two. lobes which are separated by the neutral 
sheet current boundary; the upper or northward lobe has its magnetic field 
pointing roughly toward the earth, whereas the southward lobe field points 
away from the earth. The moon is immersed in the tail about four days of each 
orbit; the moon can pass through either or both lobes. 

Global eddy current induction with the Moon in the geomagnetic tail. The 
induced poloidal field is considered to be in a vacuum in the deep-lobe regions 
of the tall. 

Transient response data, measured when moon was in the solar wind, (a) 
Nightside transient response data, showing decay characteristics of the 
radial component of the total surface field B after arrival of a step transient 

AX 

which reduces the external magnetic field radial component by an amount a B_ 

Ex 

here normalized to one. The shape of the curve illustrates time characteristics 
of the decay of the induced poloidal eddy-current field (from Dyal and Parkin, 
1973). (b) Daytime transient response data, showing decay characteristics of 
tangential components (B Ay z ) of the total surface field after arrival of a step 
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transient which increases the external magnetic field tangential component 


Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


by an amount &B , here normalized to one. Shape of the curve again 
Ji»y f z 

illustrates decay characteristics of the induced poloidal field. The overshoot 
maximum is amplified to ~ 5 by solar wind dayside compression; the 
theoretical overshoot maximum is 1. 5 for an unconfined poloidal field. (From 


Dyal et al. , 1973). 

Transient event in the geomagnetic tail, measured simultaneously by the 
Apollo 12 LSM and the Explorer 35 Ames magnetometer. Data are expressed 
in the surface coordinate system which has its origin at the Apollo 12 
magnetometer site; x is directed radially outward from the surface, while y 
and z are tangent to the surface, directed eastward and northward, respectively. 
Due to poloidal field induction in the moon, the Apollo 12 radial (x) component 
is "damped” relative to the Explorer 35 radial component, whereas the Apollo 12 
tangential (y and z) field components are "amplified 11 relative to Explorer 35 
data. Apollo 12 and Explorer 35 component data scales differ due to the 
existence of a 38 + 3 gamma remanent field at the Apollo 12 surface site. 
Electrical conductivity analysis for a transient event in the geomagnetic tail. 
Shown are data from the radial component of the event of Figure 6. Response 
to the Explorer 35 external field radial component is computed numerically for 
the conductivity profile shown in the insert and compared to the measured 
Apollo 12 response field. In the preliminary results of geomagnetic tail 
conductivity analysis, this selected conductivity profile, though not unique, 
yields a satisfactory fit of input and response data for this and twenty other 
tail events processed to date. 

Electrical conductivity profiles for the lunar interior calculated from 
measurements in different regions of the lunar orbit around the earth. The 
shaded region (3) is from nightside transient analysis in the solar wind (Dyal 
and Parkin, 1973), curve 1 is from dayside transient analysis in the solar 
wind, and curve 2 is from transient analysis in the geomagnetic tail (Dyal 
et al. , 1974). 


Temperature profiles of the lunar interior calculated from the conductivity 
profile in Figure 8. These calculations use the laboratory data for olivine 
from Duba ei al. (1972). 
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Figure 10. Magnetization induction in the moon. When the moon is immersed in a 

uniform external field H (in this case the steady geomagnetic tail field), 

a dipolar magnetization field is induced in permeable material in the 

lunar interior, with the dipole axis of_M aligned along the direction of H. 

The total magnetic field near the moon is B_ = H_+ 4 ttM . The magnetic 

permeabilities of the two layers are u and u^, and for regions outside 

the moon, u= u =1 (free space). H is measured by the lunar orbiting 
o 

Explorer 35, whereas E is measured by an Apollo Umar surface magnetometer 
(LSM). Measurements of B and H allow construction of a B-H hysteresis 
curve for the sphere, from which permeability and iron abundance can be 
calculated. 

Figure 11. Hysteresis curve for the moon. Data points are 2703 simultaneous 2-minute 

averages of radial components of the external geomagnetic field data H 

(measured by the lunar orbiting Explorer 35 Ames magnetometer) and total 

magnetic induction B = uH (B is measured by the Apollo 12 lunar surface 

magnetometer). Data points are selected from four lunations of measurements 

made when the moon was immersed in the uniform geomagnetic tail field. In 
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this low-external-field regime (~10 gammas or. 10 Oe) , the hysteresis curve 
is linear and is fitted by a least-squares line of slope 1. 008 + 0. 004. This 
slope corresponds to a whole- moon magnetic permeability of 1. 012 + 0. 006 
(from Parkin et al. , 1974a, b). 

Figure 12. Iron abundance calculated for the moon from the measured global permeability. 

The percentage of free iron is calculated for a Curie isotherm located at a depth 
of 260 km in the moon, and the total iron abundance is calculated for two lunar 
models, (1) a moon composed of free iron and orthopyroxene arid (2) a mcon 
composed of free iron and olivine (from Parkin et al. , 1974a, b). 

Figure 13. Lunar remanent magnetic fields measured at the surface Apollo 16 Descartes 
sites, (a) Photograph showing the Apollo 16 landing site, the location of the 
surface magnetometer (LSM), and the traverse positions where the portable 
magnetometer was deployed, (b) Magnitude and orientation of the measured 
vector remanent magnetic fields. 



Figure 14. Magnetic energy density versus plasma energy density at two Apollo sites 
which have different remanent magnetic fields. The magnetic energy 
density is computed from the difference between the compressed and uncom- 
pressed remanent field at the Apollo sites. Plasma energy density data 
are calculated from Apollo solar wind spectrometer (SWS) measurements. 

N is the proton number density, m is the proton mass, and v is the plasma 
bulk speed. Apollo 12 magnetometer data are plotted versus Apollo 12 
solar wind spectrometer data, while Apollo 1G magnetometer data are plotted 
versus Apollo 15 SWS data, SWS data are courtesy of C. W. Synder and 
D. R. Clay of the Jet Propulsion Laboratory. Uncompressed remanent 
field magnitudes are 38 v at Apollo 12 and 235 v at Apollo 16 LSM sites. 
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